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TWODIMENSIONAL  VISCOUS  FLOW  OF  A  COMPRESSIBLE  FLUID 
1 .  PRELIMINARIES 


The  equations  of  motion  for  a  twodimensional  flow  of  a  compressible 
viscous  fluid  are 


,  ,  13  ,3u  3v  ,32u  .  32u 

p(uux+vuy)+px  =  +  +  IfrZ  +  3^  ' 


(Navier-Stokes) 


(1.1) 


/  w  1  3  ,3u  ^  3v.  ^  ,32v  a  32v. 

p(uvx+vvy)+py  =  3^-  +  +M<^r  +  ay?) 


(pu)x 


+  (pv)  =  0  , 

y 


(continuity) 


Where  u  and  v  are  the  velocity  componants,  p  is  the  pressure,  p  is  the 
density  and  p  is  the  constant  coefficient  of  viscosity.  In  order  to 
make  the  system  determinate,  one  assumes  that  either  the  density  function 
p(x,y)  is  known,  or  that  the  equation  of  state  in  the  form  of  a  pressure- 
density  relation  is  known. 

We  shall  develop  equivalent  systems  of  equations  which  replace 
(1.1)  similar  to  those  developed  by  Martin  { 1 }  for  viscous  incompressible 
fluids,  by  Chanda  {2}  for  Cosserat  fluids  and  Nath  {3}  for  magneto¬ 
hydrodynamic  flows.  The  new  systems  of  equations  will  have  the  vorticity  w, 
speed  q,  velocity  gradient  e,  defined  as, 

(1.2) 

u>  =  v  -u  q  =  /u2+v2 ,  e  =  u  +v 

x  y,  n  x  y, 

the  pressure  p  and  density  p  as  dependent  variables.  A  system  of  curvi¬ 
linear  coordinates  defined  by  the  streamfunction  ^  and  an  arbitrary  function 
form  the  independent  variables  of  the  system.  Once  solutions  are  obtained 


for  these  equations,  the  flow  may  be  mapped  onto  the  physical  plane. 


2 


After  establishing  the  new  systems  of  equations,  we  shall  consider 
two  special  cases  of  flows,  one  with  streamlines  enveloping  an  arbitrary 
curve  and  the  second  case  with  streamlines  as  involutes  of  an  arbitrary 
curve.  In  both  cases  we  shall  show  that  the  curve  must  reduce  to  a 
single  point. 

In  order  to  arrive  at  the  systems  of  equations  mentioned  above,  we 
need  the  following  preliminary  results: 

We  first  re-write  the  system  (1.1)  as 


pqqx+px  '  ■  y  -  x 


-pva)  =  -pu)  _+4pe, 
3 


(Navier-Stokes) 


pqq  +p  +puw  =  pw  +4pe 

V  y  x  y  y 


pe+p  u+p  v  =  0 
x  y 


(1.3) 


It)  =  v  -U 

x  y. 


e  =  u  +v 
x  y 


(continuity) 

*»• 

(Vorticity,  Gradient) 


The  stream  function  i(j(x,y)  is  defined  by  the  continuity  equation  of  (1.1), 
*x  =  ~PV'  *y  =  Pu  (1.4) 

Let  <j>(x,y)  be  an  arbitrary  function  so  that  the  two  systems  of  curves 
<|>  =  const.,  i|i  =  const,  determine  a  curvilinear  coordinate  system.  We  may 
now  treat  the  variables  x,  y,  u,  v,  u,  q  and  z  as  functions  depending  on  <t> ,  ip. 
Once  these  functions  are  known  in  the  ($,i|))  -  plane,  we  may  transform  them  to 
the  physical  plane  provided  the  Jacobian 


wv*, 


in  non-zero  and  finite.  Then 


x  =  J4»  , 

<t>  y 


y  =  Jip  ,  x  = 

<t>  x  ^  y' 


y,  =  j«t  . 

\l>  x 


(1.5) 


(1.6) 
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The  first  fundamental  form  for  the  (x,y) -plane  is  given  by 
ds2  =  Ed<J>2  +  2Fd<i>dii'  +  Gdiji2 , 
where 

E  -  p  -  <(yw  G  -  and  »  - 


The  Gaussian  curvature  given  by 


K  =  i2/ r  2  _  w 

17  UE  11’ *  E 


r  2. 

12  i 


where 


r  2  __ 

li  - 


-FE.+2EF  -EE, 

<p  <p  to 

2W2 


r  2 

12 


EG.-FE, 

2W7 


(1.7) 


(1.8) 


(1.8) 


(1.9) 


r  2  = 

22 


EG , -2FF , +FG . 
- 71 - 


2W<: 


are  the  christaffel  symbols  {4},  should  vanish. 

If  6  is  the  angle  of  inclination  of  the  tangent  to  the  streamline 
=  const,  in  the  sense  of  increasing  <j>,  then 


6  -  /£  ( rifa4>  +  Tila^) , 


(1,10) 


arid 

j  (5 

Z  =  x+iy  =  /e _  {Ed$  +  (F+iJ)  dip} 

/e~ 


«,11) 


with  J  =  +W. 


The  identities 


2  E 
11  '  *2 


(1.12) 


F  P  2  E  -  2 


<«*>♦  ’  5*  riJ  -  iP  r22 


and 


<£> 

W  4> 


-  (\) ,  -  |  hi  -  £  hi  +  |  hi 

W  i|)  W  W  V} 


will  be  of  use  in  later  calculations. 


EQUATIONS 

With  the  introduction  of  vorticity  ,o,  speed  q  and  qradient  e  as 
given  in  (1.2),  the  Navier-Stokes  equations  (1.1)  may  be  written  with 
curvilinear  coordinates  i{)  as 

pqq*  =  +  %  j  •  §  - V 


4  t  G  F 

pqq.  =  — pe.  +  pw.  —  -  pw, - p,, 

>  3  iJj  <p  j  ij)  j  ip 


or,  equivalently. 


G(pq^  +  %  +  f^v  - F(pqq„, +  %  +  u  +  fv  = 


F(pqq^>  +  P<j>  +  3pEij))  +  E  (pqq,  +  P,  +  to  +  |pe  )  =  pJto  . 

(p  J  (P 

We  may  use  theorem  4.1  of  {1}  to  write  the  equation  of  continuity  as 

q  =  /e 

^  •  <2-3> 

Further  we  may  use  (1.6)  together  with  (1.4)  to  eliminate  the 

velocity  componants  in  the  equation  of  continuity  of  (1.3)  to  obtain 

e  =  Iqlq  (2.4) 

j'p  \  ’ 

v 


or,  using  (2.3) 


Thus  the  equation  of  continuity  and  the  equation  for  the  velocity 


gradient  e  may  be  replaced  by  any  two  of  the  three  equations  (2,3)  -  (2,5) 
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From  the  pair  of  equations  (1.12)  together  with  (2.3),  we  have 


(pq)  ,  =  Fill  Eri? 


(pq)  =  F  F1 2 ~ 1 '  ‘  2 2 
<l>  - 


(2.6) 


y'E  W  /E  W 

If,  using  (2.6)  and  (1.10),  the  derivatives  of  q  and  6  are  eliminated 
from  the  relation  ({l},  eg.  5.2), 


/eVIjJ  =  Fq  -E< 


we  get 


r~  EG 
/EWw  =  —  r, , 
11 


2EF_ 

p/iw 


2  .  E  2 

r  +  -  r  ^ 

12  p/m  22 


/E  , 

- 7  (Fp  -Ep  ) 

pMf  '  1  l/ 


(2.7) 


Using  (1.13),  this  equation  may  be  written  as 


W  w  *  (pw*$  “  (pW^  ' 


(2.8) 


This  equation  would  then  replace  the  vorticity  equation  of  (1.3).  we 
now  have  a  system  of  equations  which  may  be  used  to  replace  (1.3)  and  equiva¬ 
lently,  (1.1). 

Theorem  2.1.  If  the  streamlines  i|i  =  const,  of  the  flow  of  a  viscous  com¬ 
pressible  fluid  are  taken  as  one  set  of  coordinates  in  a  curvilinear  coordi¬ 
nate  system  <pr  ip  in  the  physical  plane,  the  system  (1.3)  of  five  partial 
differential  equations  for  u,  v,  w,  p  and  e  may  be  replaced  by  the  system 


G(^V%+iue^  "  F(pqVp^+U)+ =  • 

(Navier-Stokes) 

'F(pqq  *+v  +E<pqvp**»+ fvV  = 

/e  ] 

Pq  =  ^  (~)  ^  =  +EW,  (continuity,  gradient) 
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;l,{  (  -  (  7:)  ,  }, 

W  pV!  p  p  W  I,1; 


and 


K  =  W  {  (I  ‘I^V  =o. 


(vorticity) 


(Gauss) 


of  six  equations  for  E,  F,  G,  q,  u>,  P,  p  and  e  as  functions  of  <}>  and  p . 
Once  a  solution  is  found  for  this  system,  the  flow  in  the  physical  plane 
is  obtained  by  the  two  relations  (1.10)  and  (1.11). 

The  Gauss  equation  above  insures  that  the  first  fundamental  form 
given  in  (1.7)  with  E,  F  and  G  as  functions  of  (j>  and  ^  is  for  a  plane. 

We  now  proceed  to  a  second  formulation  of  the  equations  given  in 
Theorem  2.1. 

If  we  differentiate  the  first  equation  of  (2.1)  with  respect  to  p 
and  the  second  with  respect  to  4>  and  subtract,  we  arrive  at  the  equation 


I  ^  ±  q 


3 (p,q) 


♦  ) 


=  o 


(2.9) 


where  A2  is  Beltrami's  differential  parameter  of  second  order  {4}, 


.  _  ]  r.Gf.-Ff,,  .-Ff  +Ef  ,  , 

A2f  =  ^  1  (  <t>  40  +  (  $  ) 

w  <j>  W  V 


(2.10) 


and  3(p,q)  is  the  Jacobian  p,q,-p,q,  . 

-rwttt  *  *  *  * 


Using  the  equation  (2.3)  for  q,  we  may  eliminate  the  derivatives 

of  q  from  (2.9)  to  obtain 

.  )  3(1/P,E/W2) 

PWA2o>  r  ±  -  - wrff-  =  (2.11) 

We  will  now  derive  an  equation  for  the  pressure  p.  Using  (2.3),  we 
obtain  the  derivatives  of  q  and  with  the  help  of  (1.12),  (1.13)  and  (2.7) 
write  them  in  the  form. 


pqq.  =  ± 

<P 


2 

I’ll 


k 


F  1  rC 

pqq  >  +  —  1  1 

<!< 


.  ' ,  r  2 ) 

W  '  p  '  W '  ‘  1 1  W  1  7  ' 


Then  equation  (2.2)  may  be  written  as. 


-i-  -FP,+EP,  j  4  -Fe  +Ec  .  2 

pw ,  =  +  <1<  ±  — p  <J>  *  -1  I . ; 

*  W  3  W  pw  11  ' 


(2. 12) 


Gp^-Fp,  4  Ge  -Fe 
pw  =  +  4>  <P  +  JV  4>  1 1‘  ± 


pW 


r  2 

12 


We  now  apply  the  integrability  condition  w 
equation  to  obtain 


w  and  use  the  Gauss 

U’4> 


/  4  %  2 

(p  +  rpe )  =  - 
3  p 


h? 


r  2  - 

12 


(r  2)2 

12 


V]z 


i  r  '■ 
(-).  12 
P 


(2.13) 


Using  (2.8),  the  vorticity  w  may  be  eliminated  from  (2.11)  and 
using  (2.4),  e  may  be  eliminated  from  (2.13).  Then  these  two  equations, 
together  with  the  Gauss  equation  form  an  underdetermined  system  of  three 
equations  for  E,  F,  G,  p  and  p.  Once  a  solution  is  obtained  for  this 
system,  the  flow  in  the  physical  plane  is  obtained  from  equations  (1.10)  and 
(1.11).  We  now  state  the  following  theorem: 

Theorem  2.2.  The  twodimensional  flow  of  a  viscous  compressible  fluid  may 
be  described  by  first  solving  the  underdetermined  system  of  three  partial 
differential  equations  given  by  (2.11)  and  (2.13)  from  which  w  and  e  are 
eliminated  and  the  Gauss  equation  K  =  O  for  E,  F,  G,  p  and  p  and  then 
obtaining  the  angle  of  inclination  a(<p,4/)  from  (1.10)  and  x(<)>,i|)),  y ( , iJj ) 
from  (1.11)  which  will  then  implicitly  define  the  streamfunction  4>(x,y). 

We  will  now  develop  a  third  system  of  equations  for  q,  w,  p  and  p 


as  functions  of  <f>  and  We  set 


9 


E  =  vi2w  ?  +  (pqq.+P.-  ~VC  )2  , 

<}■  <?  9  3  <!> 

a.  9  4  4 

F  =  yzu)  u .  +  (pqq.+P,-  r^ic  )  (pqq  +  p,  +w-  -ye  )  , 

❖  V  $  <P  3  4>  3  ip  (2.14) 


% 

G 


w2“%+<pqVV“"  a'V 


-v  4  4 

J  =  yu  (pqq  +p  +w-  qyc) -yw  (pqq+p,-  rye  )  . 
o  9  3  y  u)  iji  A  3  9 


The  equations  (2 

.2)  may  then 

be  written  as 

E  _  F 

G  J 

W 

^  'X/ 

E  F 

'X,  'V 

G  J 

'X, 

W 

(2.15) 

with  W  =  ±J 

according  as 

W  =  ±J. 

From  (2.3) 

and  (2.15)  we  have 

W  = 

(2.16) 

p2q2W 

so  that  equations  (2.15)  yield 

Ez 

VO 

EF 

0/b 

EG 

E  =  p2q2W2  ' 

F  —  o 

'  G  =  p2q2W2  ' 

(2.17) 

which  shows  that 

if  u),  q,  e. 

p  and  p  are  known  as 

functions  of  <j>  and  ip, 

the  quantities  E,  F  and  G  may  be  calculated  from  (2.17)  which  would  then 
give  the  flow  in  the  physical  plane  from  equations  (1.10)  and  (1.11). 

The  equation  (2.5)  together  with  (2.15)  will  yield  us  the  first 


equation  of  the  new  system. 


pi  ? 

e 


{2. IS) 


The  vortieity  equation  (2.8),  which  can  now  bo  written  in  the  form 


% 

(-i)  -  (.=) 

pJ  V  pJ  V 


(0  E 

pV  j 


(2.19) 


The  Gauss  equation  may  be  obtained  from  the  integrability  condition 


n  .  ,  =  a  .  , 

where  n.  and  a  are  obtained  from  (1.10)  as 

<P  <P 


(2.20) 


J  ..  2 

%  E  ll* 


u  y 

a ,  =  —  1  f 
ip  E  12  . 


(2,21) 


To  obtain  expressions  for  the  Christoffell  symbols  and  r^2  ,  we  use 

(2.4)  to  eliminate  e  from  (2.12)  and  use  these  equations  to  eliminate  the 
derivatives  of  q  from  (2.2).  We  then  have  " 


■k  r“  =  »*“♦ +  v  K>  l  -  (v  K>  i ' 


1  r  2  ,  .  4  ,G  /IS  4  F  ,1,  1 

pw  p  iji  Vl4>  3^  V  w  <j>  r  fw  p  4>  w 


(2  v22) 


The  integrability  condition  (2.20)  then  yields  the  Gauss  equation 


in  the  form 


,  1  4  ,F  1  4  .  E  ■, 

(~2  (PA  +  rve,)- —  7772  (P,+  gut  i )—  ) 

pq  <J>  3  <J>  ^  pq*  ifi  3  ip 


ip  3  iji  a,  \fi 
J 


_{-L(P  +  Le  >£  _  L_{ p  +  *  )L  +  _£_}  .  Ip  a HUo,g) 

pq2  P((i  3^  pq^  <|>  3^  \p' ^  pq2  pq3  9(<t>,<p) 


(2.23) 


•V  - 
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Theorem  2.3.  When  the  streamliner  ’y  =  const,  of  a  steady  viscous 
compressible  flow  are  taken  as  one  set  of  coordinates  in  a  curvilinear 
coordinate  system  <p,  ip  in  the  physical  plane,  the  system  (1.1)  for  the 
velocity  componants  u,  v,  pressure  p  and  density  p  as  functions  of  x  and 
y  may  be  replaced  by  a  system  of  three  equations  (2.18),  (2.19)  and 
(2.23)  for  the  vorticity  to,  speed  a,  velocity  gradient  c,  pressure  p 
and  density  p  as  functions  of  0  and  ip.  Once  the  solutions  for  this 
system  are  found,  the  flow  in  the  physical  plane  may  be  obtained  by 
calculating  the  inclination  a(<j>,i|i)  from  (1.10),  and  y(<|>,tjj)  from 

(1.11)  with  the  help  of  (2.17). 

The  systems  of  equations  developed  so  far  can  be  made  determinate 
if  the  equation  of  state  is  unknown  and  by  a  proper  choice  of  the 
coordinate  curves  <t>  =  const.  For  example,  if  the  curves  <J>  =  const, 
are  chosen  as  orthogonal  to  the  streamlines,  the  system  given  in 
Theorem  2.3  becomes 


'V 

F  =  o  , 


% 

U)  E 
pq  2  , 


(2.24) 


{_  + 

pq  <P 


4  E, 


-  {-t-<pa+ 
pq 


-the 
3  < 


'b 

)-  +  -T-}  .  = 

j  -v  pqz  <(, 


3p  3 (m,q) 
pq7  3(4>,4') 


One  may  also  choose  <j>  =  p  to  obtain  another  system  of  equations  for 
ti),  q,  e  and  p  as  functions  of  p  and  ip  provided  the  isobars  and 


streamlines  do  not  coincide. 
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3.  Special  Cases. 

Consider  the  orthogonal  system  of  curves  consisting  of  tangent 
lines  and  the  involuter  of  an  arbitrary  curve  C.  The  system  may  be 
represented  by  two  parameters  n  =  const.,  5  =  const,  where  n  is  the 
angle  of  elevation  of  the  tangent  lines.  If  a  =  o(n)  is  the  arc- 
length  of  C  from  a  fixed  reference  point,  the  radius  of  curvature  is 
given  by 


do  =  R 
dn 

and  the  first  fundamental  form  is  given  by 
ds2  =  d?2  +  (5~o)2dn2. 


(3.1) 


(3.2) 


We  shall  first  study  flows  for  which  the  streamlines  are  tangential 
to  the  curve  C.  We  then  choose  the  second  coordinate  curves  <j>  =  const, 
as  the  involutes  of  C  so  that 


<)>  =  <H5),  =  *(n). 


(3.3) 


Comparing  (3.2)  with  (1.7)  we  get  expressions  for  t,  F,  G  and  J  as 


.  e  =  r-)2,  F  =  O,  G  =  (il sj‘‘ 


where  *  =  d 
d? 


and 


d 

dn  • 


j  =  5-a 


(3.4) 


From  (2.3),  (2.4)  and  (2.8)  we  have  expressions  for  q,  w  and  s  as 


q  = 


_ Sll 


_  -  -  ±1 

P(5-a)  ,  C-o  V?  '  - C-o  lp75-o)/n. 

The  Gauss  equation  of  Theorem  2.2  is  identically  satisfied  and  the 
equation  (2.11)  for  w  becomes 


(I) 


1  ,  ^  ■ 
w  =  -  r —  {-  v 


(3.5) 


al  (5~o) 5  +  a2(5-o)4  +  a3(5-o)3  +  *4(5-0) 2  +  a5(5-o)  +  a6  =  o, 


where 


(3.6) 
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a>  •  *  Kv  • 


=  3VJ-, 


2 

+  3(1^"  (“)  c-pi|<  'a  ’  (— )  rr  +  2*  V(i),  +  (r) 


pen 


PC 


P  cc 


p  c 


'P  Cn  , 


a,  =  5a*  V  (i)t  -  W(i)„  -  3a*'"  (i)  n  -  3»+-(i>nn  -  rt  •£>„„„ 


-4pi[>"(->  +2<K?o‘  (-)r  -  (~)  -  2v"«"(-)  , 

p  p  4  p  n  p 


(3.7) 


2  1 

a  =  (9pt|>'  ct'  -  6ptJ/"a"  -  3^*  a'  )  — 


(3,7) 


-(12pi//"0'  +  4pi> ’a")  (i)  -  6yiJj  'a  '  (— )  . 

p  n  p  nn  ' 


l  z  l 

a5  =  -15^P  ’o  ‘  (— )  ^  -  (15i|)"o  '  +  10t<>  ‘a  'a")  (~) , 


a6  =  -15^’o'  . 


Equation  (3.6)  is  an  identity  which  must  hold  everywhere,  and  it 
must  also  be  valid  on  the  curve  C,  given  parametrically  as  £  =  o(n). 

Thus  a6  =  o  which  could  happen  only  if  0’  =  0.  By  (3.1)  This  means 
that  the  radius  of  curvature  vanishes  and  the  curve  C  must  reduce  to  a 
point.  We  thus  have  the  theorom, 

Theorem  3.1.  If  the  streamlines  in  the  two  dimensional  flow  of  a 
viscous  compressible  flow  are  straight  lines,  they  must  be  either 


parallel  or  concurrent. 
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As  a  second  example  we  choose  the  involutes  of  the  curve  C  as 
streamlines  and  the  tangent  lines  as  the  orthogonal  coordinate  curves 
<*>  =  const,  so  that 

<f>  =  <Mn)  >  ^  =  iM?,)  .  (3.8) 


We  then  have  as  before 

„  (C-a)2  ___  „  1  T  _  £-a 

E  = *r  ,  F— 0  .  G  —  —  .  J  —  — '"7 

—  '  '  *  <*>  .  • 

y  y 


(3.9) 


and 


pq  =  i  ,  e  =  “-(-)  ,  w  =  -  (-)  -  ~ - r 

C-o  p  n  p  £  p(£-o) 


(3.10) 


Using  (3.9)  and  (3.10),  equation  (2.11)  may  now  be  written  as 
ai(£-o)5  +  a2  (£-o)4  +  a3  (^-o) 3  +  a4(£-c)2  +  a5(£-o)  +  a6  =  0, 
where 

^  *  3?,i,{  ♦  3i,i>ts  *  ♦(!)  )  . 


(3.11) 


»2  -  2U  4  ♦  4„hi>{  +  2p*(i)u  -  2W(i)n  42(i)E„ 


a  =  (;-)  +  yiC  (-)  +  (-)  r  -  i  (-) 

J  p  p  £  pnn  p  £nn  p  n 


(3.12) 


a,. 


=  pi(i(i)  +  P<i< (~) „  +  p^rJn0  +  pij'Aj.a'  - 


fr2c  ' 


p  no 


p'£n 


•  r  ^ 

aS  =  p{3i^(~)^a  ’  +  ^  a"}  ,  a6  =  3p  -  a  ’ 


P  0  p 


Again,  (3.11)  should  hold  on  the  curve  C,  leading  to  the 
conclusion  that  o’  =  o  and  the  curve  must  reduce  to  a  point.  Thus 


we  have 


Theorem  3.2.  The  streamlines  in  the  two  dimensional  flow  of  a 


compressible  viscous  fluid  can  be  involutes  of  a  curve  C  only  if  C 
reduces  to  a  point  and  the  streamlines  will  then  be  circles  concentric 
at  this  point. 

To  obtain  the  flow  in  the  physical  plane,  one  would  need  the 
equation  of  state  p=p(p),  Then  equation  (2.13)  together  with  either 
(3.6)  or  (3.11)  would  give  a  determinate  system  of  equations  for  the 
pressure  p  and  streamfunction  ip.  Once  these  functions  are  solved,  the 
flow  in  the  physical  plane  can  be  obtained  by  using  (3.4)  and  (3.9) 
in  (1.10)  and  (1.11) . 


Ur. 
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